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We consider the Einstein-Maxwell-dilaton system with an arbitrary kinetic gauge function and 
a dilaton potential. A family of analytic solutions is obtained by the potential reconstruction 
method. We then study its holographic dual QCD model. The kinetic gauge function can be fixed 
by requesting the linear Regge spectrum of mesons. We calculate the free energy to obtain the 
phase diagram of the holographic QCD model and interpret our result as the heavy quarks system 
by comparing the recent lattice QCD simulation. We finally obtain the equations of state in our 
model. 
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I. INTRODUCTION 

To study phase structure of QCD is a challenging and important task. It is well known that QCD is in the 
confinement and chiral symmetry breaking phase for the low temperature and small chemical potential, while 
it is in the deconfinement and chiral symmetry restored phase for the high temperature and large chemical 
potential. Thus it is widely believed that there exists a phase transition between these two phases. To obtain 
the phase transition line in the T ~ n phase diagram is a rather diflScult task because the QCD coupling 
constant becomes very strong near the phase transition region and the conventional perturbative method does 
not work well. Moreover, with the nonzero physical quark masses presented, part of the phase transition line 
will weaken to a crossover for a range of temperature and chemical potential that makes the phase structure 
of QCD more complicated. Locating the critical point where the phase transition converts to a crossover is 
another challenging work. For a long time, lattice QCD is the only method to attack these problems. Although 
lattice QCD works well for zero density, it encounters the sign problem when considering finite density, i.e. 
H 0. However, the most interesting region in the QCD phase diagram is at finite density. The most concerned 
subjects, such as heavy-ion collisions and compact stars in astrophysics, are all related to QCD at finite density. 
Recently, lattice QCD has developed some techniques to solve the sign problem, such as reweighting method, 
imaginary chemical potential method and the method of expansion in fJ-/T. Nevertheless, these techniques are 
only able to deal with the cases of small chemical potentials and quickly lost control for the larger chemical 
potential. See [l| for a review of the current status of lattice QCD. 

On the other hand, using the idea of AdS/CFT duality from string theory, one is able to study QCD 
in the strongly coupled region by studying its weakly coupled dual gravitational theory, i.e. holographic 
QCD. The models which are directly constructed from string theory are called the top-down models. The 
most popular top-down models are D3-D7 0-0 model and D4-D8 (Sakai-Sugimoto) model [1,0|- In these 
top-down holographic QCD models, confinement and chiral symmetry phase transitions in QCD have been 
addressed and been translated into geometric transformations in the dual gravity theory. Meson spectrums 
and decay constants have also been calculated and compared with the experimental data with surprisingly 
consistency. Although the top-down QCD models describe many important properties in realistic QCD, the 
meson spectrums obtained from those models can not realize the linear Regge trajectories. To solve this 
problem, another type of holographic models have been developed, i.e. bottom-up models, such as the hard 
wall model Q and the later refined soft-wall model @. In the original soft-wall model, the IR correction 
of the dilaton field was put by hand to obtain the linear Regge behavior of the meson spectrum. However, 
since the fields configuration is put by hand, it does not satisfy the equations of motion. To get a fields 
configuration which is both consistent with the equation of motions and realizes the linear Regge trajectory, 
dynamical soft- wall models were constructed by introduce a dilaton potential [lol . [Tl| consistently. Later on, 
the Einstein-dilaton and Einstein-Maxwell-dilaton models have been widely studied numerically 12-16j. By 
potential reconstruction method, analytic solutions can be obtained in the Einstein-dilaton model |13] and 
similarly in the Einstein-Maxwell-dilaton model [l^, [l^ . 

In this paper, we consider the Einstein-Maxwell-dilaton system with an arbitrary kinetic gauge function and 
a dilaton potential. A family of analytic solutions are obtained by the potential reconstruction method. We 
then study its holographic dual QCD model. The kinetic gauge function can be fixed by requesting the meson 
spectrums satisfy the linear Regge trajectories. We study the thermodynamics of the Einstein-Maxwell-dilaton 
background and calculate the free energy to obtain the phase diagram of the holographic QCD model. By 
comparing our result with the recent lattice QCD simulations, we interpret our system as the model for the 
heavy quarks system. We finally compute the different equation of states in our model and discuss their 
behaviors. The behavior of the equations of state is consistent with our interpretation of the heavy quarks. 

The paper is organized as follows. In section II, we consider the Einstein-Maxwell-dilaton system with a 
dilaton potential as well as a gauge kinetic function. By potential reconstruction method, we obtain a family 
of analytic solutions with arbitrary gauge kinetic function and warped factor. We then fix the gauge kinetic 
function by requesting the meson spectrums to realize the linear Regge trajectories. By choosing a proper 
warped factor, we obtain the final form of our analytic solution. In section HI, we study the thermodynamics of 
our gravitational background and compute the free energy to get the phase diagram. From the phase diagram, 
we argue that our background is dual to QCD with heavy quarks and interpret the black hole phase transition 
as the deconfinement phase transition in QCD. We finally plot the equations of state in our background to 
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compare with that in QCD. We conclude our result in section IV. 



II. EINSTEIN-MAXWELL-DILATON SYSTEM 



We consider a 5-dimensional Einstein-Maxwell-dilaton system with probe matters. The action of the system 
have two parts, the background part and the matter part, 



S = Sb + Sr, 



(2.1) 



In string frame, the background action includes a gravity field , a Maxwell field and a neutral dilatonic 
scalar field 6^, 



Sb 



1 



~g e 



-'2<t's 



Rs 



fs 



(2.2) 



where fs {<j>s) is the gauge kinetic function associated to the Maxwell field and Vg {(j>s) is the potential of 
the dilaton field. One should note that the function fs {(j>s) is a positive-definite function. The explicit forms 
of the gauge kinetic function fs {(t>s) and the dilaton potential Vs {4>s) are not given ad hoc and will be solved 
consistently with the background. 

The matter action includes the flavor fields (^Aj^,Aj}), which we will treat as probe, describing the degree 
of freedom of mesons on the 4d boundary. 



S„ 



1 



fs 



{Fl + 



R) ' 



(2.3) 



where G5 is the coupling constant for the flavor field strength F^ij= d^A^ — d^A^^ and /s (<^s) is the gauge 
kinetic function associated to the flavor fleld. The gauge kinetic function of the flavor field in the matter action 
is not necessary to be the same as that of the Maxwell field in the background action |2.2| in general, but here 
we set them equal for simplicity. 

In the above, we have defined our model in string frame in which it is natural to write the boundary conditions 
as we will see when solving the background in the next section. However, to study the thermodynamics of 
QCD at finite temperature, it is convenient to solve the equations of motion and study the equations of state in 
Einstein frame. To transform the action from string frame to Einstein frame, we make the following "standard" 
transformations, 



V 8 ^^"^ 
The background and the matter actions become, in Einstein frame, 

1 



, /,(^,)=/(^)eV3-, Vsm^e V3V(0) 
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'9 



'9- 



:/(</>) 



+ Ft, 



(2.4) 

(2.5) 
(2.6) 



where we have written the fiavor fields A^ and A^ in terms of the vector meson and pseudovector meson fields 
V and V, 



A^ ^V + V, A^^V-V. 



(2.7) 
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The equations of motion can be derived from the actions (|2.5I) and (|2.6I) as 



^ d(f) Adcf) 
= 0, 



R 



0, 
0, 

/(0) 



2 43m-^' + {^V',^^v} 



1 2 

df_,(j)di,(j) - -gf,i, {d(t>) - g^^uV 



(2.8) 

(2.9) 
(2.10) 

(2.11) 
(2.12) 



In the next section, we will solve the above equations of motion under some physical boundary conditions and 
constraints. 



A. The gravitational Background 



In this section, we will solve the background of the Einstein-Maxwell-dilaton system defined in the last 
section. We first turn off the probe flavor field V and V in the equations of motion (|2.8II2.12|) . which reduce to 



. ^ dV df 



d(j) 4 df 

V,, [/(<^)^^^n = 0, 



1 2 

d^(j)dy(j) - -g^„ {d(j)) - gf,„V 



(2.13) 
(2.14) 
(2.15) 



Because we are interested in the black hole solutions, we consider the following form of the metric in Einstein 
frame, 



-g{z)di 



> (z) , ^ At (z) 



(2.16) 
(2.17) 



where z = corresponds to the conformal boundary of the 5d spacetime. We will set the radial L of AdS's 
space to be unit in the following of this paper. 



Using the ansatz of the metric, the Maxwell field and the dilaton field (|2.16[ [^T7|) . the equations of motion 
and constraints for the background fields become 



'-+3A'-^ 



25 



z'^g 







A' 

z 6 



3g' 6\ ^, 1 fig' 4\ , g" , e^^V 



A" + 3A" + 



z \ 2g z j 6g 3z^g 



0, 



= 0, 



(SA'- J) g'-e-2^zVAf = 0, 



= 0. 



(2.18) 

(2.19) 

(2.20) 
(2.21) 

(2.22) 



We should notice that only four of the above five equations are independent. In the following, we will solve 
the equations (12. 19112. 22l) . and leave the equation (|2.18l) as a constraint for a consistent check. 
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To solve the background, we need to specify the boundary conditions. Near the horizon z = z^f , we require 

At{zn)=g{zH)^Q, (2.23) 

due to the physical requirement that A^A^^ — q'^^AqAq must be finite at z = z^f . 

Near the boundary z — > 0, we require the metric in string frame to be asymptotic to AdS^^ thus 

1 



dsl_ 



^ (z 0) dx'^dx" = — [-dt'^ + dz^ + dx^] . 



(2.24) 



This boundary condition, in Einstein frame, becomes, 



^(0) 



0(0), .9(0) = 1. 



(2.25) 



Since we do not assume the form of the dilaton potential V {(j))^ which should be solved consistently from 
the equations of motion, we will treat the dilaton potential as a function of z, i.e. V {z), when we solve 
the equations of motion. With the above boundary conditions (|2.23l) and p.25p . the equations of motion 
()2. 19112.22]) can be analytically solved as 



</)'(z) 
At (z) 

9iz) 



'-6{ A" - A'^ + -A' 



re ^"^dyj^^ ^^^-^"^ ^^■'^ 



dy, 



dy n —r-:dx 



irv'^-'^dyJl-^dx 



V{z) = -iz^ge-^^ 



A" + 3A'- 



25 



A' 



2<7 



65 



(2.26) 
(2.27) 

(2.28) 
(2.29) 



where we have used the boundary conditions to fix most of the integration constants. The only undetermined 
integration constant yg will be related to the chemical potential fi in the following way. We expand the field 
At (z) near the boundary at z to get 



At (0) 



-1 



\j;"y'e-^^dyjl^^dx 



y .1 



1 



6-4(0)/ (0) 

Using the AdS / CFT dictionary, we can define the chemical potential in our system as 



(2.30) 



-1 



'^dyJl^^dxJo 



dy, 



(2.31) 



in which yg can be solved in term of the chemical potential ^ once the gauge kinetic function / (z) and the 
warped factor A{z) are given. Put the solution ([2.26112.2^ into the constraint (|2.18p . it is straightforward to 
verify that the above solutions are consistent with the constraint. 

We note that the solutions (|2. 26112. 29| depend on two arbitrary functions, i.e. the gauge kinetic function 
/ (z) and the warped factor A (z). Different choices of the functions / (z) and A (z) will give different physically 
allowed backgrounds. Thus we have just found a family of analytic solutions for the Einstein-Maxwell-dilaton 
system. We will use the freedom of choosing functions / (z) and A (z) to satisfy some extra important physical 
constraints. 
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B. Vector Meson Spectrum 



In a theory with hnear confinement like QCD, the spectrum of the squared mass of mesons is expected 
to grow as n at zero temperature and zero density. This is known as the hnear Regge trajectories [19|. In 
the method of AdS/QCD duahty, this issue was first addressed in Q by modifying the dilaton field in the IR 
region with a term, i.e. the soft- wall model. In the term was added by hand to the dilaton field. 
It means that the fields configuration used in the soft-wall model is not a solution of the Einstein equations. 
Dynamically genera ting the z'^ term by consistently solving the Einstein equations has been considered in 
several later works jT7r[20| by including a proper dilaton potential. At finite temperature and density, the 
temperature dependent meson spectrum has been studied in [2T| - [26| with the AdS thermal gas background 
replaced by the charged AdS black hole background. 

In the previous section, we consistently solved the equations of motion (|2. 18112.2^ for the Einstein-Maxwell- 
dilaton system. The analytic solutions depend on two arbitrary functions, the gauge kinetic function / (z) and 
the warped factor A (z). In this section, we will study the meson spectrum in our background and constrain 
the functions of / (z) and A (z) by requesting the vector meson spectrums satisfy the linear Regge trajectories 
at zero temperature and zero density. 

We consider a 5d probe vector field V whose action has been written down in (|2.6p . 

5™ = / d'xV^^F^ (2.32) 



167rG5 „ 

To get the meson spectrum, we study the vector field V in the charged AdS black hole background which we 
have obtained in the previous section, 



2Aiz 



'g{z)df + —— + dx^ 
9{z) 



(2.33) 



= At (z) dt. (2.34) 

The equation of motion of the vector field V has been derived in (|2.10p as 

[/ {4>) i^r] = 0. (2.35) 

Following |9!], we first use the gauge invariance to fix the gauge Vz — 0, then the equation of motion of the 
transverse vector field [d^V^i = 0) in the background (|2.33l) reduces to 

ivV, + V'^ + ( ^ + 4 + ^' - - ) X^/ = 0' (2-36) 
5 \g f zj 

where the prime is the derivative of z. We next perform the Fourier transformation for the vector field Vi as 

F,(x,z)= /-^e*-i;.(z), (2.37) 
J (27r) 

where k — (j-u,p) and the functions Vi (z) satisfy the eigen-equations 
Redefining the functions Vi (z) with 

1/2 

= Xij,, (2.39) 



brings the equation of motion (j2.38l) into the form of the Schrodinger equation 

-.2 



g 



< + C/(z)V^, = ( (2.40) 
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where the potential function is 



U(z) 



2X2 



(2.41) 



In the case of zero temperature and zero chemical potential, we expect that the discrete spectrum of the vector 
mesons obeys the linear Regge trajectories. At = 0, the metric of the zero temperature background (thermal 
gas) coincides with the black hole metric in the limit of zero size, i.e. zjj — >■ cx), which corresponds to g (z) = 1. 
In the zero size black hole limit, the Schrodinger equation reduces to 

+ C/ (z) V', = m2^„ (2.42) 

where — — = —up- To produce the discrete mass spectrum with the linear Regge trajectories, the 

potential U (z) should be in certain forms. Following [§], a simple choice is to fix the gauge kinetic function 
as / (z) = e^'^^ -^(2)^ then the potential becomes 



2 2 

c z . 



The Schrodinger equations (|2.40p with the above potential (j2.43p have the discrete eigenvalues 



4cn, 



(2.43) 



(2.44) 



which is linear in the energy level n as we expect for the vector spectrum at zero temperature and zero density. 
At finite temperature and finite density, 5 (z) 7^ 1, the masses of the vector mesons solved from the Eq. (I2.40p 
will depend on the temperature and density. For the case of small enough temperature and density, Eq. (|2.40p 
can be solved perturbatively to get the mass shift from the linear Regge trajectories HJ-flll- For large 
temperature and density, the method of spectral functions is useful. The study of temperature and density 
dependent vector mass spectrum is in progress. 

Once we fixed the gauge kinetic function / = e^'^^ -MA ^ the Eq. (|2.31l) can be solved to get the integration 
constant yg in term of the chemical potential /i explicitly as 



(1 -e^^^y 



/q " y^e -^^dy 2cy? J~" y^e ^^dy ^ ^ 

Put the integration constant yg back into the solution (|2. 26112. 29|) . we can finally write down our solution as 



(A'(z) 

At [z) 
9{z) 
V{z) 



' -&[ A" - A'^ + -A' 



Jo 



-3z2ge-2^ 



~^^dy 
A" + 



2c^2 



1 - 

3A'2 



/J 

J Zf 



2h 3 - 

ye 



^^dy /; 



y^e 



V _ 6\ , _ 1 / V _ 4 
2g z) z \ 2g z 



9_ 

6g 



dy 
dy 



y^e' 



3A 



dy 



(2.46) 
(2.47) 
(2.48) 
(2.49) 



Now we have fixed all the integration constants by either satisfying the boundary conditions (|2.23p and (I2.25P 
or relating to the chemical potential /x. The final solution (12.46112.4^1) depends only on the warped factor A (z). 
The choice of A{z) is arbitrary provided it satisfies the boundary condition (j2.25l) . In the next sections, we 
will make a simple choice of A (z) and use it to study the phase structure in its holographic QCD model. 



III. PHASE STRUCTURE 



We will study the phase structure for the black hole background which we obtained in the last section 
(|2. 46112. 4^ . The phase transitions in the black hole background correspond to the phase transitions in its 
holographic QCD theory by AdS/QCD duality. 
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A. Fixing the Warped Factor 



To be concrete, we fix the warped factor A (z) in our solution in a simple form as 

A{z)^~^z'-bz\ (3.1) 

where the parameters b and c will be determined by later. It is easy to show that this choice of A (z) satisfies 
the boundary condition (|2.25p by Eq. (j2.46p . There are many more complicated choices for the function of 
A (z), but we will show that our simple choice already educes abundant phenomena in QCD. 

Once we choose the function of A (z) up to the parameters b and c, our solution (|2. 46112. 49|) is completely 
fixed. Expanding the dilaton field and the dilaton potential near the boundary z — 0, respectively, 



0(z) = 2v^z+— ^ ^^z^ + ---, (3.2) 

9v3c 

V = -12-18cz^ + ■■■ , (3.3) 

we can write the dilaton potential in terms of the dilaton field cj) as the expected form from the AdS/CFT 
dictionary, 

y = -12 + + - ,A^3. (3.4) 

The conformal dimension A = 3 satisfies the BF bound 2 < A < 4 implying that our gravitational background 
is stable. Furthermore, the dilaton satisfying the BF bound corresponds to a local, gauge invariant operator 
in 4d QCD possibly. One should note that the dilaton in our work is different from that in the references 
p^| - [l8| . where the dictionary between the dilaton and gauge invariant operator in field theory are not clear. 
In our case, we combine the effects of quarks and gluon which absorbed by dilaton potential. 

We fix the parameter c by fitting our mass formula^ = 4c (n + 1) to the lowest two quarkonium states^, 

mj/^ = 3.096Gey, m^, = 3.685Gey, (3.5) 



For c ~ lAGGeV^, we have 



mi = SMGGeV, to2 = 3.731Gey, (3.6) 



which are consistent with the experimental data within 1%. The parameter b will be determined later by 
comparing the the lattice QCD result of the phase transition temperature at zero chemical potential. 



B. Black Hole Thermodynamics 



Using the black hole metric we obtained 

2A{z 



-g{z)dt^ 



dz^ 
9iz) 



dx^ 



(3.7) 



^ We shift our mass formula 1 12. 441 1 by 1 to make n = 1 correspond to the lowest quarkonium states, J/V"- 

^ According to the analysis of the phase structure of our background that we will see later in this section, we would like to 
interpret our holographic QCD model to describe the heavy quarks system. 
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it is easy to calculate the Hawking-Bekenstein entropy 




0.6 0.8 1.0 1.2 

(a) (b) 

FIG. 1: The temperature v.s. horizon at the different chemical potentials ^ — 0, 0.35, 0.5, 0.714, 0.9Gey are plotted. 
We enlarge the rectangle region in (a) into (b) to see the detailed structure. At ^ = 0, there is a global minimum Tmin5 
for < jj, < jj,c — 0.714:GeV , the minimum becomes local and eventually disappears for /i > fic- 



We plot the temperature T v.s. horizon Zh at different chemical potentials in FIG.[T] At ^ = 0, the temper- 
ature has a global For Zh > ^min, the black hole solutions are thermodynamically 
unstable. Below the temperature Tmin, there is no black hole solution and we expect a Hawking-Page phase 
transition happens at a temperature Thp ^ Tmin where the black hole background transits to a thermal gas 
background. For < ^ < fic — O.TlAGeV , the temperature has a local minimum/maximum temperature 
I)i mm/7)i max at rH = Tmin/rmax and dccrcascs to zero at a finite size of horizon. The black holes between 
rmin and rmax are thermodynamically unstable. We expect a similar Hawking-Page phase transition happens 
at a temperature Thp > 7)imin- In addition, since the thermodynamically stable black hole solutions exist 
even when the temperature below T^min, we also expect a black hole to black hole phase transition happens 
at a temperature Tbb G [J)iminj 'Tjimax]: where the large black hole transits to a small black hole. The values 
of Thp and TpB will determine the true vacuum state, thermal gas or small black hole, in which the system 
will stay eventually. Finally, for /i > /Xc, the temperature monotonously decreases to zero and there is no 
black hole to black hole phase transition anymore'^, which implies that there is a second order phase transition 
happens at ii — ^c, i.e. the critical point. 

To determine the phase transition temperatures Thp and Tbb, we compute the free energy from the first 
law of thermodynamics in grand canonical ensemble, 

F ^ e~Ts- fip. (3.10) 

Changes in the free energy of a system with constant volume are given by 

dF = -sdT - pdfi. (3.11) 



^ There could still be a Hawking-Page phase transition at some temperature for the case of fi > fic, but we will show later that 
the black hole solution is always thermodynamically favored in the case. 
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At fixed values of tfie chemical potential fi, the free energy can be evaluated by the integral 27 1 



F = - 



/ 



sdT. 



(3.12) 



We can fix the integration constant in the above integral (|3.12p by considering the zero chemical potential case. 
At /z = 0, the metric of the zero temperature background (thermal gas) coincides with the black hole metric 
in the limit of zero size, i.e. zj^ od, where we expect that the free energy of the black hole background also 
coincides with the free energy of the thermal gas background which we can choose to be zero. Thus we require 
F (z// — >• oo) = and obtain that 



FIG. 2: (a) The free energy v.s. horizon at the different chemical potentials /i = 0,0.35,0.5,0.71,0.9. At /i = 0, we 
normalized that the free energy vanishes when zh oo. For ^ < fic — 0.71, the free energy has a maximum; for 
/i > jj,c, the free energy becomes monotonous, (b) The free energy v.s. temperature at the different chemical potentials. 
At ^ = 0, we normalized that the free energy vanishes when zh — > oo. For fj, < fic, the free energy has a maximum; for 
fj, > fj,c, the free energy becomes monotonous. 

With the choice of A (z) in (|3.ip . the integral in (|3.13p can be performed to get the free energy of the black 
hole at fixed chemical potentials. We plot the free energy v.s. horizon at different chemical potentials in FIG. 
[2] We have normalized the free energy of the black holes by requiring it to vanish (or equal to the free energy 
of the thermal gas) at zjj — > oo. Therefore, the black hole background is favored for the negative value of the 
free energy and the thermal gas background is favored for the positive value of the free energy. At ^ = 0, the 
free energy starts from a large negative value at small zh to a positive maximum value and then decreases 
to zero at zjj — > oo. The free energy intersecting the x-axis implies that there exists a Hawking-Page phase 
transition from the black hole to the thermal gas background at Thp- For < fi < fic, besides a maximum 
value, the free energy has a local minimum value, which implies a black hole to black hole phase transition. 
For fi > fJ-c, the free energy becomes monotonous and no phase transition exists. 

The phase structure is more transparent in the plot of free energy v.s. temperature in FIG. O At fi — 0, 
the free energy increase from a negative value with a large temperature to zero at T = Thp where the black 
hole transits to the thermal gas background which is thermodynamically stable for T < Thp- At finite /i, the 
free energy behaves as the expected swallow-tailed shape. For < fi < fic, the curve of free energy intersects 
with itself at T = Tbb where the large black hole transits to the small black hole background. We found that 
the free energy of the black hole is always less than that of the thermal gas, i.e. Fbiack hole < -^thermal gas = 0. 
Therefore the thermodynamic system will always favor the small black hole background other than the thermal 
gas background when T < Tbb- When we increase the chemical potential fi from zero to fj,c, the loop of the 




(3.13) 




(a) 



(b) 
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FIG. 3: The phase diagram in T and /i plane. At small p, the system undergoes a first order phase transition at 
finite T. The first order phase transition stops at the critical point (ncTc) — (0.714Gel/, 0.528GeV), where the phase 
transition becomes second order. For ji > Hc, the system weaken to a sharp but smooth crossover. 



swallow-tailed shape shrinks to disappear at ^ = /ic, where the background undergoes a second phase transition. 
For fi > fic, the curve of the free energy increases smoothly from higher temperature to lower temperature. 

We plot the phase diagram of our holographic QCD model in FIG. |3l At /i = 0, the system undergoes a 
black hole to thermal gas first order phase transition at T = T^p. For < fi < n^., the system undergoes a 
large black hole to small black hole first order phase transition at Tbb- The phase transition temperature Tbb 
approaches to Thp at — > that makes the phase diagram continuous at /i = 0. By comparing the phase 
transition temperature at /i = to the lattice QCD simulation of Thp — Q.GGeV in ;28;], we fix the parameter 
b ~ 0.273GeV^. The first order phase transition stops at the critical point (ncTc), where the phase transition 
becomes second order. For fi > fic, the system has a sharp but smooth crossover. 

The phase diagram of our holographic QCD model in FIG. [3] is not expected from the common picture 
of QCD phase diagram, in which the transition is crossover for the small chemical potential /i and sharpens 
to the first order phase transition beyond a critical fic- Then question is how to interpret our result? To 
explain our result, we need to look at the phase structure of QCD more carefully. Lattice QCD provides many 
useful information about the phase diagram of QCD at least at /i = 0. For finite density, lattice QCD suffers 
the well-known sign problem. However, several perturbative methods such as reweighting, complex chemical 
potential and expansions in /x/T etc. have been developed in lattice QCD to compute the physical quantities 
at finite density. Despite that these methods only work for small chemical potential /i, they can help us to 
understand the whole picture of QCD phase diagram. 

FIG. m shows the schematic phase transition behavior of three flavor QCD for different choices of quark 
masses at ^ = [l|. We can see that, the order of the phase transition at ^ = depends on the quark masses'*. 
There are two limit regions where the phase transition is first order at /i = 0. One of them is near the chiral 
limit, i.e. mu = rud = rus = 0. The other is near the decoupling limit, i.e. rriu — rrid — vtis — ^ oo. Thus there 
are two possibilities to interpret the phase diagram of our holographic QCD model. 

The first possibility is to interpret our result as for the light quarks (near the point of chiral limit) and the 
phase transition is the chiral symmetry breaking phase transition. For a not very large chemical potential 
/X, lattice calculation shows that the chiral critical line separating the crossover from the first order phase 
transition regions expands as increasing the chemical potential /i to form a chiral critical surface. The sign 
of the curvature of the critical surface is crucial to determine the phase structure at finite fi. In FIG. [5l the 
chiral critical surfaces in the case of positive and negative curvature are showed. In the interpretation of the 



* Most of the lattice results prefer tfiat the physical mass point locates in the crossover region, but this is not completely confirmed 
yet j9j. 
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FIG. 4: Schematic phase transition behavior of three flavor QCD for different choices of quark masses at zero density 
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FIG. 5: The chiral critical surfaces in the case of positive (left) and negative (right) curvature from lattice QGD [lt|. 



light quarks, the mass point locates near the origin. The first order transition will be preserved for any finite 
fj, if the curvature of the chiral critical surface is positive, while it will becomes a crossover at a finite chemical 
potential if the curvature of the chiral critical surface is negative. Recent lattice calculations shows that the 
curvature of the chiral critical surface is more like negative [11, that is consistent with the phase diagram 
of our holographic QCD model. 

The other possibility is to interpret our result as for the heavy quarks (near the point of decoupling limit) and 
the phase transition is the deconfinement phase transition. Similar to the chiral critical line, lattice calculation 
shows that the deconfinement critical line separating the crossover from the first order phase transition regions 
expands as increasing the chemical potential fi to form a deconfinement critical surface. FIG. IHl shows the 
recent result of the deconfinement critical surface from lattice QCD [28l |. In our interpretation of heavy quarks, 
the mass point locates near the infinite mass point. The first order phase transition at ^ = will becomes a 
crossover at a finite chemical potential fic and it is consistent with the phase diagram of our holographic QCD 
model. The phase transition for heavy quarks has been extensively studied recently by using lattice QCD 
method [ElliiO-lil. 

Among the two possible interpretations of our holographic QCD model, we will focus on the heavy quark 
interpretation in the rest of our paper. The reason is that, for the light quarks, the lattice QCD simulation 
has almost confirmed that the physical point of light quarks mass locates in the crossover region at zero 
chemical potential in FIG. ID While for the heavy quarks, the lattice QCD does not give any constraint yet 
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FIG. 6: The deconfinement critical surface from lattice QCD [28 



and the physical point of the heavy quarks mass has great possibility to locate in the first order region at zero 
chemical potential. Furthermore, if we re-examine the gravity side more carefully, we realize that we did not 
take into account the backreaction from the matter fields when we solved our gravitational background. We 
only consider the baryon number chemical potential, but not the dynamical quarks. This means that we are 
considering the quenched limit of the heavy quarks. 

To interpret our result as for the heavy quarks with the deconfinement phase transition, there is still a 
problem in the gravity side. It is widely believed that the deconfinement phase transition in the field theory 
side is dual to the Hawking-Page phase transition in the gravity side. Hawking-Page phase transition is the 
transition between black hole and thermal gas backgrounds. However, in our gravity background, the phase 
transition is between a large black hole and a small black hole backgrounds for non-zero chemical potential. 
Thus it sounds inconsistent to consider that the black hole to black hole phase transition in the gravity side 
is dual to the deconfinement phase transition in QCD. Our resolution for this problem is that, although it is 
thermodynamically stable, the small black hole is dynamical unstable. On the other hand, the gauge group in 
realistic QCD is SU [N] ^ SU (3). Thus we have to consider the finite iV effect in the gravity side. If we still 
require that the gravity limit held, this is merely to consider the finite string coupling constant by including 
the loop correction, i.e. the quantum effect. It is well-known that small black holes are quantum mechanically 
unstable and will quickly evaporate away by quantum radiation. Therefore, right after the phase transition 
of a large black hole to a small black hole, the small black hole will continue to evaporate away quickly to a 
thermal gas background. In this sense, the black hole to black hole phase transition can be interpreted as the 
deconfinement phase transition in the dual QCD theory. 

In our heavy quarks interpretation, an important question is where the critical point is, i.e. the point 
where the first order phase transition cease to become a crossover in the phase diagram FIG. [3] Locating 
the critical point is a crucial job to understand the phase structure of QCD. As we have seen in FIG. [U 
the critical point is the point where the self-intersection disappears. We thus obtain the critical point as 
(/ic,T'c) — (0.714Gey, 0.528G'eV"). To justify the critical point we got, we expand the the function At {z) in 
Eg. (12.471) near the boundary z — >• as 



At iz) 



1 



from which the baryon p density can be read off as 



2c/i 



2c/i 

1 - e'^^H 



(3.14) 



(3.15) 



We plot the baryon density p v.s. the chemical potential p, in (a) of FIG. [T] For T < Tc, the baryon 
density p is single- valued which indicates that there is no phase transition. While for T > Tc, p is multi- valued 
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(a) (b) 

FIG. 7: (a) The baryon density p v.s. the chemical potential p,. The curve is single-valued for T > Tc, while it is 
multi-valued for T < Tc. At T — Tc, the slope becomes infinite at = ^c- (b) The specific heat Cv v.s. Temperature. 
The negative value of Cv corresponds to the thermodynamically instability. For jj, < ^c, Cv is negative for a range of 
T. For fi > He the specific heat is always positive implying no phase transition. 



which indicates that there is a phase transition at certain value of the chemical potential fi. At the critical 
temperature T = the slope of the p — fJ- curve becomes infinite at the critical chemical potential pc- The 
behavior of the baryon number p (/i) near the critical temperature Tc is consistent with the result we obtained 
from the free energy. 

The susceptibility is defined as 




(3.16) 



which is just the slope in (a) of FIG. [T] For T < Tc, the susceptibility is always positive, x > implies that 
the black hole with any chemical potential value is thermodynamically stable. On the other hand, for T < Tc, 
X could be negative for a range of p where the black hole is thermodynamically unstable. At T = Tc, % — >■ oo 
at p = p^c which indicates that a phase transition happens around there. 

The similar behavior of the susceptibility x can be seen from the plot of the specific heat Cy v.s. the 
temperature T in (b) of FIG. [71 where the specific heat Cy is defined as 

We note that, in the Cv — T diagram, the negative value of the specific heat corresponds to the thermody- 
namically instability. For p > p,c, the specific heat is always positive, Cy > implies that the black hole with 
any temperature is thermodynamically stable. While ioi {) < p < pc, Cy could be negative for a range of T 
where the black hole is thermodynamically unstable. Ai p — Q, there is a minimum temperature Tmin for the 
black hole solutions where the specific heat diverges. The Hawking-Page like phase transition happens at a 
temperature slightly above Tmin at Thp- 



C. Equations of State 



The speed of sound is defined as 



n dhiT 
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FIG. 8: The squared of speed of sound v.s. temperature. For /i < fic, (?s is negative for a range of T implying 
dynamical instability and there is a phase transition. For ^ > /Xc, (?s becomes smooth and always positive. At /x = Hc, 
Cs touches the a;-axis aX T — Tc where the phase transition transfers to a crossover. 



FIG. m plots the squared of speed of sound v.s. the temperature T . For < ^ < Ug, the speed of sound 
is imaginary for a range of temperature, indicating a Gregory-Laflamme instability (35l , ISq ] . This is related to 
the general version of Gubser-Mitra conjecture [37h39|, i.e. the dynamical stability of a horizon is equivalent 
to the thermodynamic stability. In our system, the negative specific heat implies thermodynamically unstable. 
While the imaginary speed of sound implies the amplitude of the fixed momentum sound wave would increase 
exponentially with time, reflecting the dynamical instability. Roughly speaking, Cy < is equivalent to < 
in our system. For /i > /Xc, the speed of sound behaves as a sharp but smooth crossover. At the critical point 
/i = /ic, a second order phase transition happens where goes to at the critical temperature but never 
becomes negative. In all the case, approaches the comformal limit 1/3 at very high temperature as expected. 

We plot other equations of state in FIG.[9l The entropy of our black hole solution has been calculated in (|3.8p 
and is plotted in (a) of FIG. [9l At /i = 0, there is a minimum temperature Tmin for the black hole solutions. 
The black hole solutions with very low entropy and high temperature always have negative specific heat and 
are thermodynamically unstable and the black hole will transit to the thermal gas through a Hawking-Page 
phase transition. For < fi < fic, the entropy is multi- valued for a region of temperature which indicates a 
phase transition between high entropy and low entropy black holes. For /i > /ic, the entropy is single- valued 
and there is no phase transition. The similar phase behaviors have been discussed in [12] for a holographic 
QCD model with different values of parameters tuned by hand. 

The pressure p and the energy e can be calculated from the free energy as 



(b) and (c) of FIG. [5] plot the pressure p and the energy e v.s. temperature T. For some particular /i and 
T, the pressure becomes negative which seems not make sense for a real physical system. However, we see 
that the pressures for all the thermodynamically favorite backgrounds are always positive. Both the pressure 
and the energy increase with the chemical potential, that pushes the phase transition temperature Tbb to 
the smaller values for growing p. Our results are consistent to the recent lattice results with finite chemical 
potential [ioj . 

We finally plot the trace anomaly e — 3p v.s T in (d) of FIG. [HI At /i = 0, the sharp peak with a infinite 
slope edge indicates a phase transition. With the growing chemical potential /i, the peak becomes wider with 
decreasing height. The phase transition turns to be weaker and weaker for increasing fi and eventually becomes 
a crossover. 



P 



F, e = F + sT. 



(3.19) 
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( c) ( d ) 

FIG. 9: (a) The entropy s v.s. temperature, (b) The pressure p v.s. temperature. (c)The energy e v.s. temperature, 
(d) The trace anomaly e — 3p v.s. temperature. 

IV. CONCLUSION 

In this paper, we studied the Einstein-Maxwell-dilaton system. We obtained a family of analytic black hole 
solutions by the potential reconstruction method. We then studied the thermodynamic properties of the black 
hole backgrounds. We computed the free energy to get the phase diagram of the black hole backgrounds. In 
its dual holographic QCD theory, we are able to realized the Regge trajectory of the vector mass spectrum 
by fixing the gauge kinetic function. We then discussed the possible interpretations of the phase structure 
that we obtained from the gravitational background by comparing the lattice QCD simulations. We argued 
that the heavy quarks interpretation is more favored in our system. We calculated the equations of state in 
our holographic QCD model. We found that our dynamical model captures many properties in the realistic 
QCD. The most remarkable feature of our model is that, by changing the chemical potential, we are able 
to see the conversion from the phase transition to a crossover dynamically. We identified the critical point 
in our holographic QCD model and calculated its value. As the authors knowledge, our model is the first 
holographic QCD model which could both dynamically describe the transformation from the phase transition 
to the crossover by changing the chemical potential and realize the linear Regge trajectory for the meson 
spectrum. 

Since we interpret our holographic QCD model as a heavy quarks system, it would be interesting to perform 
a lattice simulation on the equations of state of a heavy quarks system and compare with our results. On the 
other hamd, there are many future directions one can study. For example, the most interesting issue is to find 
an appropriate warped factor such that one can obtain a phase diagram similar to the common QCD phase 
diagram. Another interesting issue is to incorporate the chiral symmetry breaking by introducing a scalar 
coupled to the flavor fields and take the backreactions of flavor fields into account. It is also interesting to 
compute the linear quark-antiquark potential and expectation value of Polykov loop. One can compute the 
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meson spectrum and determine the quarkonium dissociation temperature in our background. One can also 
compute the various transport coefficients hke shear visocisty, bulk viscosity and so on. It is also interesting 
to compute the critical exponents of various physical quantities near the critical point. Some of these issues 
are in progress. 
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